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in 

I The inverse proximity effect, i.e. the induction of a magnetic moment in the superconductor 

'•—J ■ in superconductor-ferromagnet (S/F) junctions is studied theoretically. We present a microscopic 

approach which combines a model Hamiltonian with elements of the well established quasiclassical 
theory. With its help we study systems with arbitrary degree of disorder, interface transparency and 
thickness of the layers. In the diffusive limit we recover the result of previous works: the direction 
of the induced magnetization M is opposite to the one of the F layer. However, we show that in the 
ballistic case the sign of M may be positive or negative depending on the quality of the interface 
and thickness of the layers. We show that, regardless of its sign, the penetration length of the 
magnetic moment into the superconductor is of the order of the superconductor coherence length, 
which demonstrates that the effect has a superconducting origin. 
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I. INTRODUCTION 



The proximity effect in multilayered structures consisting of superconductors and normal metals has been studied 
intensively in the last decades. Since the sixties it is well known that superconducting correlations can penetrate 
into a normal region changing the physical properties of the latter 0. More recently special interest has been paid 
I to the study of the proximity effect in superconductor-ferromagnet structures. The presence of an exchange field 
. acting on the spin of the electrons leads to interesting new physics, as for example oscillations in the superconducting 
O ' critical temperature and in the density of states (DoS) of the ferromagnet as a functions of its thickness 
, The interplay between superconductivity and magnetism also leads to a change of the sign of the Josephson current 

in S/F/S structures 5, 6] (for a review see Q)- On the other hand, recent theoretical studies have focused on the 
inver se p roximity effect, i.e. the induction of a magnetic moment in a superconductor in contact with a ferromagnet 
[1,0,0,01. This phenomenon has been analyzed in Refs. [1,0 for diffusive S/F junctions within the quasiclassical 
] Green functions formalism. The results of Ref Q demonstrate that the direction of the induced magnetic moment is 
opposite to the one in the ferromagnet. The authors of Ref.Q also shown that this effect is related to the induction 

j of a triplet component (with the zero projection S'^ = of the total spin) in the pairing amplitude. 

The interest in the inverse proximity effect has been increased with the development of experimental techniques like 
ly-j neutron reflectometry 0| and muon spin rotation which allow to determine accurately the spatial distribution 
of magnetic moments. For example experiments on multilayered system consisting of the high Tc superconductor 
YBa2Cu307 and the ferromagnet L2/3Ci/3Mn03 layers have shown an induced magnetic moment in the S layers with 
a sign opposite to the magnetization of the ferromagnet |l2j| . However, a conclusive explanation of this experiment 
Q ■ has not been given yet. 

More recently the inverse proximity effect was studied theoretically in atomic-scale F/S/F trilayers 0|. It was 
shown that if the magnetization of the F layers points in the same direction, the S layer acquires a spin polarization 
parallel to the one in the F layers. At first glance this result contradicts the one of Ref.j^], but until now no 
detailed comparison between these two works has been done. For completeness we mention that in Ref.0 the induced 
magnetization in a clean S/F bilayer was computed for highly polarized ferromagnets and a leakage of the magnetic 
moment into the superconductor over atomic distances was obtained. In Ref.0| the induced magnetic moment in a 
I clean superconductor in contact with a magnetic wall was studied. 

^ In this work we present a microscopic model which allows to study the inverse proximity effect for systems with 
an arbitrary degree of disorder, interface transparency and layer thickness. The model is based on the Hamiltonian 
approach of Refs. [TsIIi^IitI . Within this model we show how disorder can be introduced by an appropriate averaging 
procedure in such a way as to recover the quasiclassical theory in the diffusive limit. 

The main results of the paper can be summarized as follows: For an atomic-size bilayer in the pure ballistic limit, 
in which the component of the momentum parallel to the plane of the layers is conserved, the induced magnetic 
moment has a positive sign in accordance with the results of Ref. jlO| | . However this is only valid for low values of the 
interface transmission coefficient r. When r is of the order of unity the sign of the induced magnetization may be 
negative. These results can be understood in terms of the behavior of the bilayer Andreev states as discussed in Sect. 
III. When disorder is included and the interface becomes diffusive we recover the result obtained in Ref. , according 
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to which the induced magnetic moment in S is negative regardless of the interface transparency. We also consider 
the intermediate situation and show that one can go smoothly from the ballistic to the diffusive regime. In section 
IV we consider the case of a S layer with thickness much larger than the superconducting coherence length ^g. We 
show that the sign of the induced magnetic moment in the superconductor depends on the thickness of the F layers 
as well as on the transmission of the interface. We also show that the penetration of the magnetic moment into the 
superconductor is in all cases of the order of ^s- This demonstrates that the effect is governed by superconducting 
correlations and is not just a "magnetization leakage" over atomic distances into the S region. 

The rest of the paper is organized as follows: in Sect. II we introduce the model Hamiltonian. Details of the model 
are presented in the appendix, where we also show the equivalence between our model and the quasiclassical approach 
in the appropriate limit. In section III we study the inverse proximity effect in a S/F bilayer of atomic-size, while in 
section IV we consider layers of arbitrary thickness. Finally section V is devoted to the concluding remarks. 

II. THEORETICAL APPROACH 

In order to describe a S/F layered structure we use a lattice model Hamiltonian as in Ref. . This can be written 
as H — Hs + Hp + Hs f where 

Hs ^ ^ es»a+^a»,^ + ^ A,(a+^a+_ +aj_aj,+) + ^ tsijiaf^^aj^^ + a+^a-j^a) (1) 

2,(7 — it i <2,J>(7 — ± 

± i <'t,j>.cr— lb 

are the terms describing the uncoupled superconducting and ferromagnetic regions. Indexes i,j label different sites 
along the system. The coefficients tij are the hopping amplitudes between neighboring sites. The site energies are 
measured from the Fermi level. is the local order parameter and hi is the local exchange field in the ferromagnet. 
Throughout this paper we assume that A (h) is a given parameter which is constant and finite in the superconductor 
(ferromagnet) and is equal to zero in the ferromagnet (superconductor). 

Finally Hsf is the term which describes the coupling between the superconducting and ferromagnetic regions, 
which is given by 

Hsf = ^■S-P^^J (aj>ai:<T + a^^aj,^) . (3) 

i&S,j&F,a=± 

Written in this way the Hamiltonian can be used to describe many different physical situations. In particular we 
would like to consider not only ballistic but also diffusive systems. One can include disorder either in the diagonal 
terms, i.e. introducing a random correction Sci for the site energies, or Stij for the hopping terms. One can also model 
a disordered interface by adding a random term Stij to the hopping elements of Eq. 

The Green functions Gn{k), which are 4x4 matrices in the particle-hole (Nambu) (8) spin space should be determined 
from the Dyson equation (see for example Refs. [TtLITsI ') 

Gs(F) = 9SiF) + gsiF){k)^s{F)Gs{F) , (4) 

where gs(F) is the Green function at the interface on the S(F) side for the uncoupled regions. As discussed in the 
appendix^ the self-energy t^s(F) in the ballistic case is given by 

^S{F)if^) = tsFTi9F(S){k)Ts. , (5) 

while in the case of a diffusive interface has the form 

^iF)=5t^r^{Gp(^S)±)Tz, (6) 

where < ... > denotes integration over parallel momentum k. The problem is then reduced to a closed set of self- 
consistent equations which can be solved as discussed in the next sections. 

With the help of this model we can also describe a diffusive metal by assuming that also the hopping terms between 
the planes inside the S and F layers are randomly distributed. As we show in the appendix after averaging over 
disorder our model leads to a discretized version of the well known Usadel equation ^19j . Thus, our model reproduces, 
in the diffusive limit, the results of Rcfs.|3.l20| for the induced magnetization. 



FIG. 1: Geometry of the S/F structure. 



An intermediate situation between the purely ballistic and the diffusive limits can be analyzed by introducing a 
parameter a controlling the proportion of the hopping term which is random, i.e. Ib — odsp and 5t — [1 ~ a)tsF 
in such a way that a = 1 (a = 0) corresponds to the ballistic (diffusive) case. In a general intermediate case the 
self-energy can be written as a sum of a ballistic S^^^^ and a diffusive S^^^,^ contribution, calculated according to 
Eqs. 10) and ®. 



III. THIN S/F BILAYER 



Let us consider two thin S/F layers, such that the unperturbed Green functions are given by Eq. HA1|I . Our aim is 
to calculate the magnetic moment M induced in the superconductor, which is given by 

M^^JiBTY,TI[{Gs+{k,uJ)-Gs-{k,u)) . 
As it was mentioned above we assume a linear dispersion relation (flat bands) . This means that 




where v is the normal density of states at the Fermi level. Notice that in the wide-band approximation that we are 
using here the magnetization arises from changes in the density of states close to the Fermi level. 

The exact normal Green functions Gs,f± for spin up and spin down can be written in the more general case as (see 
Eq.©) 



Gs± = 
Gf± — 



9s± 



9f± 



^F± ^ ^F± 



(7) 
(8) 



The self-energies S''''^ correspond to the terms defined in Eas. (|5l6|l . We thus obtain the following set of coupled 
equations 



< <-'S± > = ~^\,4 , 0+2 , o . _ , n , ^2 I A2 , ,2\> (9) 



h% + 2tl{LjA± -elT ASt^ < Fs± >) + Dsil^i + + St^ < Fs± >)' 



<Gf±> = -^( , , ^2^....' ". ,\^.2 , ,2 , ..4.^ .^ ) (11) 



. p . ^ I Tt%A + SePs < Fs± > 

^t% + 2t%{ujA±-elTASt^ <Fs±>) + Ds{Al+el + St^<Fs±>)' ^ ' 

where we have defined the following quantities: D± = oj\ -t- e^; Ds = -I- e| -I- A^; Q,± = lo + iSt^ < Gf± >', 
A± = ±A -I- 6t'^ < Ff± > and A± = uj± + i5t^ < Ff± >■ This set of equations must be solved in a self-consistent 
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way. However for the pure ballistic case, when 5t = 0, the equations in both layers are decoupled. For example it is 
easy to show that 



Gs4 



iuj{uj1 + el) + it^uj+ 

" {4~4M~4) 



(13) 



where 



ei,2 



- A2 + 2<2 ± /(cj2 - ^2 + A2)2 _ 4^2 [(^ _^ ^^)2 ^_ A2] 



(14) 



For Gs~ we obtain a similar expression substituting cj-i- by Now the integration over e can be performed by 
closing the integration path in the upper half-plane and using the residue theorem. 

In Fig. [2 we show the dependence of the induced magnetic moment on the interface transparency, for very low 
temperature (T ^ A) and different values of the parameter a and of the exchange field h. M is given in units of the 
Pauli paramagnetic moment Mp defined as 

Mp = 2iiBvh . 

The interface transparency is characterized by the transmission coefficient r , which within this model is given by 
T = 4(fsi?)2/(l -(- (tsi?)^)^- Throughout this paper all energies are given in units of ts — tp. 

According to Fig. |3 in the pure ballistic case (a = 1) the dependence M(t) has a nonmonotonic behavior: in 
accordance with the results obtained in Ref. for low transparency, the induced magnetic moment is positive. 
However, at some r (when tsp is of the order of h) M changes to negative values. If the exchange field is large enough 
this change of sign would never take place and the magnetization remains positive for all values of r. 

Figure[21also shows that by decreasing a, i.e. by increasing the disorder at the interface, the curves become flatter. 
In the limiting case of a diffusive interface, i.e. when a = 0, our model leads to the same result as the one obtained 
using the quasiclassical approach in Ref.jg]- This is not surprising since, as it is shown in the appendix, our model 
in the diffusive limit {a = 0) leads to the Usadel equation and therefore all conclusions of Ref.Q also hold here. In 
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FIG. 2: Induced magnetization M/Mp as a function of the transmission r for (a) h = 0.2 and (b) h = 0.4 and different values 
of a. We have chosen A = 0.1 and T = 10"*. All energies are given in units of ts = tp. 
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particular the physical picture given in Ref. [g of Cooper pairs sharing their electrons between both the S and the F 
layer, may explain the negative induced magnetic moment. From Eas. (|9ll2f) one can easily check that in the Usadel 
limit {a — 0) and for low transparencies the first correction to the density of states, in both layers, is proportional to 
St^, i.e. to the square of the interface transparency. This indicates that in this limit the inverse proximity effect arises 
from higher order processes requiring the induction of a non- vanishing paring amplitude in the ferromagnet. In fact, 
the Green functions in the F region have a BCS-like form with an induced minigap in the DoS We are thus 

dealing with two coupled condensates with Cooper pairs sharing the electrons between the S and the F layer. This 
coupling tends to align the spins of both regions in opposite direction. 

It would be interesting to find a physical argument for explaining the observation of a positive magnetization in the 
ballistic regime. In this case one can check from Eqs. (|9I12II that if a = 1 and for low transparency, the first correction 
to the induced magnetic moment is of order t'gp and is positive (c/. Rcf. lOj). However, when t becomes larger 
the induced magnetization reaches a maximum and then decreases down to negative values. One can understand 
this behavior in terms of Andreev states induced in the DoS. In FigOwe sketch the position of these states in the 
spectral density. For values of tsp smaller than h the state corresponding to the majority spin is located below the 
Fermi energy and thus, the induced magnetization has a positive sign. The corresponding density of states is shown 
in the inset of Fig. |2| where a peak in the spin-up channel (solid line) can be identified. By increasing the interface 
transparency and when tsF ^ h the Andreev states for both spin orientations eventually cross. A further increase of 
T leads to a negative magnetization, since now the Andreev state below the Fermi level corresponds to the spin-down 
orientation (the d-^"-^ " "--^^ '—^ t^-^ ^ 1—-^ ). 
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FIG. 3: Sketch of the position of the Andreev states in the spectral density for the majority (t) and minority (J,) spins as a 
function of the hopping parameter tsF- The three panels at the bottom show the DoS of the spin-up (solid line) and spin-down 
(dashed line) electrons for three different values of tsF- At low transparencies and below the Fermi level there is a peak in 
the spin-up channel (left panel). Increasing the value of tsF this peak moves to the right while the peak in the spin-down 
channel moves to the left. At some tsp of the order of h both peaks coincide (middle panel). Further increase of tsF leads to 
a spin-down peak below the Fermi level (right panel). 



The S/F bilayer of atomic thickness analyzed so far should not be considered as a mere idealized example but tt 
could also describe the situation in real hight Tc materials consisting of alternating magnetic and superconducting 
layers. 



IV. S/F LAYERS OF ARBITRARY THICKNESS 

Let us consider a superconductor of thickness ds ^ (where = 2afs/7rA is the coherence length) in contact 
with a ferromagnet of arbitrary thickness which is the situation in most experiments [13 . |2^ . This case can be 
analyzed theoretically considering the superconductor as a semi-infinite system. We shall first consider the ballistic 
regime in which electrons are specularly reflected (fc is conserved) at the outer surface of the ferromagnet. In Fig. 0| 
we show the magnetic moment induced in the superconductor at the S/F interface as a function of the transmission 
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FIG. 4: The induced magnetization M/Mp at the S side of the interface as a function of the transmission coefficient r for 
different values of the F-layer thickness dp- The latter is given in units of tt^f. The S layer is assumed to be semi-infinite. We 
have chosen h = 0.2, A = O.f and T=0. . 

coefficient t and for different thickness dp oi the F layer. In the tunnelling limit the effect is very weak as the induced 
magnetization is proportional to t'gp. However, in the region of high transparency an appreciable magnetic moment 
can be induced in the superconductor. As it can be observed in Fig. 0] its sign depends on the ratio between the 
thickness of the F layer and the ferromagnetic length defined as = atp/T^h: while for thin layers {dp < '^S.f) is 
negative, by increasing the F layer thickness it becomes positive reaching a maximum for dp ^ 27r^i?. For dp 3> S,p 
the induced magnetization tends to zero. Since the structure is purely ballistic we expect that interference effects take 
place leading to a nonmonotonic behavior. The change of sign can be understood within the description of Andreev 
bound states similar to the one depicted in the previous section. According to the calculations of Ref. ^3 if is 
of the order of unity, the crossing of the Andreev peaks associated to the majority and minority spins, and thus the 
change of sign of the induced magnetization, occurs at dpo ^ 4^f- This is in agreement with the present results. 

A natural question that arises at this point concerns the penetration length of the induced magnetization. In Fig. 
|31 we show the spatial dependence of M in the superconductor for two different thicknesses of the F layer and the 
same values of h and A as in the previous figure. We have chosen dp = n^p and 2tt£_p in order to show that whatever 
the sign of the induced magnetization is, the penetration length is of the order ^5. This demonstrates that the effect 
described here is due to the superconducting proximity effect and it is not just a magnetization leakage over atomic 
distances as in Ref. 

Finally we analyze the case in which the two layers are connected by a disordered interface. In Fig. we show 
the dependence of M/Mp on the transmission for different values of the parameter a and dp = 2tt£^p. We observe 
that by increasing a, the induced magnetization is reduced. However, even in the case of a fully disordered interface, 
a finite magnetization of the order of 10~^A/p for large transparency is obtained, demonstrating that the effect is 
robust against disorder. 

V. CONCLUSIONS 

We have studied the inverse proximity effect in S/F structures with arbitrary thickness, type of interfaces and 
degree of disorder. For this purpose we have presented a microscopic model which combines a model Hamiltonian 
with elements of the quasiclassical theory. We have shown that the sign and magnitude of the magnetic moment 
induced in the superconductor is highly dependent on the system parameters. In the diffusive limit, where our model 
leads to a discretized version of the well known Usadel equation, the induced magnetization in S and the one in F 
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FIG. 5: The spatial dependence of the magnetization induced in a semi-infinite S layer in contact with a F layer of thickness 
TT^F and 27r^F- We have chosen ft = 0.2 and A = 0.1 . 

have opposite directions, in agreement with Ref. 0. However, we have shown that the induced magnetization can 
exhibit either sign when the system deviates from the diffusive regime. 

We have first considered thin layers and shown that in the purely ballistic case and for low transparency of the 
interface the induced magnetization is positive. By increasing the transmission the magnetization reaches a maximum 
and then decreases to negative values. The change of sign takes place when the coupling parameter between the 
layers is of the order of the exchange field. This behavior can be understood in terms of localized Andreev states and 
their evolution when varying the coupling strength. The crossover between the ballistic and diffusive limits was also 
analyzed. 

Finally we have studied a bulk superconductor attached to a ferromagnetic layer of arbitrary thickness dp . We have 
shown that for high interface transparency and small dp the induced magnetization M is negative. For increasing 
thickness and when dp is of the order of S^p a change of sign in the induced magnetization takes place. In both cases 
(positive or negative sign of the magnetization) M penetrates into the superconductor over a distance comparable to 
the superconducting coherence length demonstrating that the effect is long-ranged. 

It is worth mentioning that band structure effects as well as a d-wave symmetry for the superconducting order 
parameter can be straightforwardly included in our Hamiltonian approach. This will be the object of future works. 

We expect that experimental techniques such as neutron rcflcctomctry or muon spin rotation used recently in Refs. 
[T^ITsI may help to confirm our predictions. 
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APPENDIX A: THE SELF-ENERGY TERM FOR DIFFERENT TYPES OF INTERFACES 

Let us consider the purely ballistic case which occurs when the mean free path of electrons is much larger than the 
characteristic dimension of the system. As we are considering layered structures, with translational invariance parallel 
to the interface {y — z plane in Fig. it is convenient to perform a change of basis from the local one used in writing 
the model Hamiltonian H1I3|) into a mixed basis of Bloch waves in the direction parallel to the interface and localized 
states in the perpendicular direction. Thus, we subdivide the layers into planes parallel to the interface, separated 
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FIG. 6: Magnetization in the superconductor as a function of r for different type of interfaces: a = 1 corresponds to the 
balhstic while a = to the diffusive one. The other values chosen for a are 0.5 and 0.75. Here dp = Ztt^f, h — 0.2, A — 0.1 
and T = 0. 

by a lattice spacing a |23|. We label these planes by integer numbers n. The new basis is given by the state vectors 
fc, n >~ J2jen e'*^^^ \ j >, where k is the momentum parallel to the plane of the junction, while Rj is the position of 
the site j within the plane. A similar model was considered in Ref. (23l] . 

After changing the basis, the field operators in Eas. (|ll3() are now labelled by k and n. We assume that the 
superconducting order parameter is homogeneous in the superconductor and therefore = A. In principle the 
dispersion relation efe may be non-linear. However in the present work we assume that all energies involved in the 
problem (A, h, etc.) are much smaller than the Fermi energy and thus we can assume a linear and isotropic dispersion 
relation, i.e. e/c ~ vpik — kp), where vp is the Fermi velocity. A generalization of this model beyond this quasiclassical 
approximation and including a more realistic band structure could be implemented straightforwardly. 

Since the system is purely ballistic, the parallel momentum fc is conserved across the interface. This means that 
the hopping matrices (Ts, Tp, Tsp) in the new basis are diagonal, i.e. 

< k'n\T\kn + 1 >= Skk'tk- 

For simplicity we shall assume that the lattice model is such that tk — J2<ij> e.'^^^^^tij do not depend on fc. This is 
consistent with the approximation of a linear dispersion relation introduced above. Thus, the self-energy term in the 
Dyson equation Eq.Q is given by 

In order to calculate the interface Green functions gs{F) of the uncoupled regions one needs to know the functions 
^S(F) uncoupled planes. These are given by [a 

(0) ( 9fiF)+ + fs{F)+ ^ \ (AU 

where the normal and anomalous components, g^g^p-^j- and Z^*^^-)^, for spin up ("+") and down("-") are known and 
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given by 

= Lo^ + ^^\el ^^^^ 
- -^^^ , (A4) 

where are the Pauh matrices in particle-hole space, ijj± = ui^ih, and w = 7rT(2n + 1) are the Matsubara frequencies. 
It is clear that the anomalous component of the unperturbed Green function equals zero in the F layer. 

The purely ballistic case is unrealistic for describing many experimental situations. Even when the uncoupled S 
and F layers are clean enough, some degree of disorder at the interface is in general unavoidable. We will show that 
the inverse proximity effect strongly depends on the type of interface and the degree of disorder. 

We start by analyzing the case of a disordered interface between the S and the F layer. In this case the layers are 
connected by random hopping terms 6tij. These random terms are assumed to follow a Gaussian locally correlated 
distribution, i.e. 



StijStim = St SimSjl- (A5) 

The observable qnt^T^^i^'^t^c cird p^lr'nlci+cirl Vnr cnrtirci rrin cr mror rlicnrrlor rrmfi mireif innc T'tiiQ cnrcircj (Tin o" T~>2'Q(2g(^-Qj'g rCStOrCS 

the translational ike in the ballistic 

case (Eq. Q) bi e disordered case 

we use the self-cc g. □ 24, 25] . The 

thick lines in this topping processes 

and the dashed t the self-energy is 
then given by 



where < ... > dei 




FIG. 7: The self-energy within the self-consistent Born approximation. The thick lines in this figure correspond to the fully 
dressed Green function, the crosses represent the hopping processes and the dashed semicircle indicates averaging over the 
parallel momentum k . 



APPENDIX B: COMPARISON WITH QUASICLASSICS 



According to our model a diffusive metal can be also described by a random coupling Stij between the layers with 
a distribution given by Eq. (|A5|) . The Green function at the plane n is determined by the Dyson equation 

(Bl) 



Gn gn + gnSt^ < TsGn+lTs > Gn + QnSt'^ < TsGn-lT^ > Gn 



where gn denotes the Green function for the uncoupled plane (cf. Eqs. (jA2IA4ll ) and < ... > the average over 
momentum. Multiplying this equation by g~^ from the left and by T3 from the right and subtracting the conjugated 
equation multiplied by g~^ from the right and by T3 from the left one obtains 



rag, 



^1 Gr, 



Gn 



, Gr, 



(B2) 



where Gn — G„T3, E„ = < G„ > and [,] denotes a commutator. 

As the dependence of g~^ on the momentum only appears in the term proportional to T3 which vanishes when 
subtracting the two Dyson equations, the integration over momentum can be done easily and is equivalent to substitute 
G by < G >. One obtains 



T-35„ \ Gn = (< Gn+1 > + < Gn-1 >) < G„ > -t^ < G„ > G„+i > < G„_i > 



(B3) 



10 



Note that by subtracting Eq. ljBl|) from its conjugate the inhomogeneous term also cancels out. Thus, we need 
a second equation in order to determine the Green functions. As in the derivation of the quasiclassical equations, 
in the wide-band approximation considered in this work, the Green functions satisfies the normalization condition 
< G >^= c, where c is a constant '2^1 . 

One can associate the terms in the r.h.s of Eq. IjBSp with the second spatial derivative of a continuous model by 
means of the identification 

< G„+i > + < G„_i > -2 < Gn > _ 2 r f'DA\ 

From the normalization condition follows that the second derivative of vanishes and therefore 
Thus, using Eqs. (|B3IB5|) we obtain following equation for g 

[Tw'\g] = -2t\^vd4gd,g) 

Identifying the factor t^a?v with the diffusion coefficient one sees that this expression is equivalent to the well known 
Usadel equation . 
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